Application of higher order holonomy corrections to perturbation theory of cosmology 
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Applying the higher order holonomy corrections to the perturbation theory of cosmology, the 
lattice power law of Loop Quantum Cosmology, ft oc p 13 , is analysed and the range of /3 is decided 
to be [-1,0] which is different from the conventional range —0.1319 > /3 > —5/2 [J]. At the same 
time, we find that there is a anomaly free condition in this theory, and we obtain this condition in 
the vector and tensor mode. We also find that the nonzero mass of gravitational wave essentially 
results from the quantum nature of Riemannian geometry of loop quantum gravity. 
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I. INTRODUCTION 



The spacetime metric of Big Bang cosmology is ho- 
mogeneous and isotropic Priedmann-Robertson- Walker 
(FRW) metric. However, this model is just an approxi- 
mation of "zero order" universe If we only focus on 
the FRW metric, we will ignore many of interesting things 
in the universe such as galaxy clusters, galaxies, stars, 
etc. So it is necessary to introduce the inhomogeneous 
and anisotropy perturbation to describe these things Q . 
On the other hand, the effects of quantum gravity should 
be significant in the very early universe. Therefore, it is 
interesting to study possible quantum gravity effects in 
cosmological perturbation theory. 

At present, the problem of finding the quantum the- 
ory of the gravitational field is still open. One of the 
most active of the current approaches is loop quantum 
gravity. Loop quantum gravity (LQG) is a mathe- 
matically well-defined, non-perturbative and background 
independent quantization of general relativity. Its cosmo- 
logical version, the loop quantum cosmology (LQC) [7] 
have achieved many successes. A major success of LQC 
is the resolution of the Big Bang sing ularity d S, E3; 
this result depends crucially on the discreteness of the 
spacetime geometry. With such a result, the big-bang 
singularity will be avoided through a big-bounce mech- 
anism in the high energy region. In addition, LQC can 
also setup suitable initial conditions for successful infla- 
tion fl3l . fl4| as well as possibly leaving an imprint in the 
cosmic microwave background [T3 |. 

In LQG, spacetime is quantized. The geometric opera- 
tors, such as the area operator and the volume operator, 
have discrete eigenvalues. So there is the smallest area 
gap A [TBI, [l6j]. In LQC, the coordinate size of a loop is 
ft, 2 , ft is the function of p — a 2 (where a is the scale factor 
of the universe.), i.e. ft = ft(p). In the early literature 
9, 10], the work always base on the simplest choice of 
ftip) = /io = const. However, this form can lead to some 
unusual features. As pointed out in [12[ that the choice 
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of ft{p) = const can lead to the Big Bounce occurs at 
classical matter density like water, so he suggest to se- 
lect the function as ft(p) oc v^ 1 ! 2 . From this time forth, 
in most of current works [ljj, °^ V 1 ^ 2 has been ap- 
plied. And it was shown that the choice ft(p) oc p -1 / 2 
is physically and mathematically consistent [l^. Up to 
now, however, there is still no theory to decide the func- 
tion of ft(p). As research continues, there may be some 
other form of ft(p) can give the better physics. There- 
fore, to find out the form of this function has theoretical 
significance. 

An ansatz for the form of this function can be taken as 
ft(p) oc pP . In [1], the range of /3 has been decided to be 
—0.1319 > j3 > —5/2. However, it is just the conclusion 
of the first order holonomy corrections. If we want a 
more accurate determination of the range of /3, we must 
consider the higher order corrections. 

Even in the case of homogeneous and isotropic models, 
the quantum equation of state is very difficult to analyze. 
Fortunately, there is a powerful tool, i.e. effective theory, 
which allows us to include loop quantum effects by correc- 
tion terms in equations of the classical type 17J. There 
are two types of quantum corrections that are expected 
from the Hamiltonian of LQG. One correction arises for 
inverse powers of the dcnsitizcd triad, which when quan- 
tized becomes an operator with zero in the discrete part 
of its spectrum thus lacking a direct inverse. The other 
comes from the fact that a loop quantization is based 
on holonomies, i.e. exponentials of the connection rather 
than direct connection components |18j . 

In LQC, there is no well-defined quantum operator cor- 
responding to c = jk. So we should find a well-defined 
operator to replace it. The conventional way is replacing 
the c by sinftc/ ft. 

The application of inverse triad corrections and con- 
ventional holonomy corrections on the scalar mode of 
perturbation can be viewed in jl9j |. the vector mode in 
|2fj| and the tensor mode in [l[. 

In this paper, we focus on the higher holonomy cor- 
rections rather than the conventional correction. We ap- 
ply these higher corrections to the vector and the tensor 
mode and see whether the mass of gravitational wave is 
the nature of discrete geometry. We will also analyses 
the range of /3 with high order holonomy corrections. 
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This paper is organized as follows. At first, the per- 
turbed variables are introduced in Sec. [TT| And then in 
Sec. IIII| we apply the high order holonomy corrections 
to obtain the effective Hamiltonian constraint. Detailed 
analysis of the range of (3 will be given in Sec. IIV1 Section 
Mis our discussion. 



II. BACKGROUND AND PERTURBED 
CONSTRAINT 

In Ashtekar's formalism of general relativity [2l|, |22| |. 
the spatial metric as a canonical field is replaced by the 
dcnsitized triad Ef , defined as 



Ef 



dot 



(4) 



and the spin connection T l a which is 



1 



r o = -e« k e) ( d [a e b] ~r - Cfc 



1 ' -eieldiA 



(1) 



(2) 



The canonical variables are densitized triad Ef and 
Ashtekar connection A l a — T^+jK^, where K % a is extrin- 
sic curvature and 7 is Barbero-Immirzi parameter. 

The canonical variables reduced to spatially homoge- 
neous and isotropic cosmology are 



E?=p8f, Ki = kSl P o =0. 



(3) 



They are background variables, and the perturbation will 
be added based on these variables. 

In the perturbation theory, we denote the background 
variables by a bar: 

E?=p5?, r a = o, K = ksi, n = Vp; N a = o, 

where p = a 2 , and the spatial metric is q a b = a 2 5^. We 
use conformal time in this paper, so we set N = a. 

The canonical variables are the perturbation densitized 
triad Ef and and Ashtekar connection A l a , which are 

Ef = pSf+SEf, Al = K +1 Kl = 7 kSi+(STi + ySKl) , 

(5) 

where 8Ef and 8K % a are small perturbation around ho- 
mogeneous variables. 

As described in [l|, Gil the symplectic structure 
splits into two parts: one for the background variables 
and the other for perturbations, i.e., 



{~k,p} = 



8ttG 



and 



{SKi(x),SE^(y)} = 8nG8 3 (x,y)S b J*. 



(6) 



(7) 



Here, G is the gravitational constant and Vq is a fiducial 
volume. 



In vector mode, the gravity part of perturbed Hamil- 
tonian constraint (up to quadratic terms) is [20] 



H G [N] 



1 



16ttG y s 



d 6 xN 



^ 8E?5E*5*SV 



2k 



(8) 



On the other hand, when we introduce the inhomoge- 
neous perturbation, the diffeomorphism constraint does 
not vanish any more. So the gravitational part of diffeo- 
morphism constraint is changed into [20( 



D G [N a ] = 



1 



8nG 



d 3 x8N c [-p(d k 8K k ) - k8 k (d d 8E$] . 

(9) 



Using Eqs.© and (J7]), we can testify the following rela- 
tion easily 



{H G ,D G } = 0. 



(10) 



Similarly, in tensor mode, the gravity part of the per- 
turbed Hamiltonian constraint (up to quadratic terms) 
is 



H G [N] 



1 



i6ttG y s 



<TxN 



2k 



8E<8Ej8 k X 

2 p3/2 



-Vp(8KlSK k 8 c k 8*) - ^(8E?SKi) 



1 



53/2 



(8 cd S^ k 8 e fd e E]dfEi) 



(11) 



where SEf — —^phf, here h l a := S lb h a b, and h a b is the 
symmetric metric perturbation field. It is transverse and 
traceless, i.e. it satisfies d a h a i, — and S ab h a b = 3. 



III. VECTOR AND TENSOR MODE WITH 
HIGHER ORDER HOLONOMY CORRECTIONS 

A. Higher order holonomy corrections 

Instead of the conventional way of introducing the 
holonomy corrections, in this article, we focus on the 
higher order holonomy corrections [23j j . 

At first, let's consider the Taylor series 



sin 1 x = ■ 



(20! 



„2;+i 



l= ^2« («r(2Z+l) 
for — 1 < x < 1 and setting x — sin (Ja^k) , we have 
(20! 



(12) 



1 k = 7,Y, 



A^2 2i (ny(2i+i) 



[sin (fijk)] 



21+1 



(13) 
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This inspires us to define a nth order holonomized con- 
nection variable as 



» 



1 " 



(20 



M^2 2 < (liy (21 + 1) 



[sin (pfyk)] 



21 + 1 



(14) 



which can be made arbitrarily close to jk as n —> oo. We 
can see that is a function of the holonomy sin [pjk) 
and the discreteness variable p. Therefore, we can replace 
jk by c h to implement the underlying structure of LQC. 

When n = 0, cl°^ = sin (pjk) //i is the same with the 
conventional holonomy corrections. 

There is an ambiguity in this replacement. If we set 
x = sin (mpjkj in Eq. (|12[) . where m is an arbitrary con- 
stant, Eq.([T3"|) changes to 



_^ oo 

nwyk = — • 



(20! 



[sin (mp'yk)] 21+1 , (15) 



M^2 2 ' («r(2/ + l) 
and we have 

(20! 



7& 



77 .//. * ' 



£jj 2 2 ' (Z!) z (2^ + 1) 



[sin (mp'yk)] 



2/+1 



(16) 



So we can define a more general nth order holonomized 

■ 0. 

h- 

(21)1 



connection variable c^: 



» ._ 1 V 



[sin (m/X7fc)] 2 ' +1 , (17) 



m/i ^ 2 2i (i!)" (2^ + 1) 



where m is an ambiguity parameter. 

The Poisson brackets between the canonical variables 
and the are 



» 



and 




8ttG 



8irG dp 
3pV dp 



cos (pjk) & n (pjk) , (18) 



i(p^k)(8 n (p-yk) k 



» 



where 



®n(Plk) =Y1 



m 

o2 2 ' (*0 ; 



[sin (/l7fe)] 



2/ 



(19) 



(20) 



From [23j we can see that the role of higher order holon- 
omy corrections is like a filter, which excludes the impact 
of human factors on the theory and leaves a pure quan- 
tum effect. 



B. Vector mode 

In classical perturbation theory of cosmology, the 
jauge invariant variables of the vector mode will decay 



quickly. Therefore, there is a little role of the vector mode 
perturbation for a universe Q. However, once we intro- 
duce the quantum correction, we must consider whether 
the perturbation theory is anomaly free 20]. The re- 
quirement of anomaly free can reduce some ambiguities 
of LQC. Inserting the higher holonomy corrections in 
Eq-©! we can obtain the effective gravity part of the 
perturbed Hamiltonian constraint 



H*[N] 



d xN 




16ttG 7e 



+Vp(6KlSK k d S c k 5f) 



2 r (n) 

_rnh SE c SK 3 

VP 7 3 



(21) 



~~ (n) 

General speaking, we should replace all the jk by c mh . 
But in order to get a homogeneous limit which agreement 
with what has been used in isotropic models, we set the 
parameter m in the first term to equal one [l|. The pa- 
rameter m in the last term should lead to an anomaly-free 
constraint algebra, so we do not fix it at first. In the fol- 
lowing discussion, we will determine the right value of m 
in the last term by requiring an anomaly-free constraint 
algebra in the presence of quantum corrections. 



In homogeneous and isotropic model, there is no dif- 
feomorphism constraint. So the algebra of constraints is 
closed. When we consider the inhomogeneous perturba- 
tion, the diffeomorphism constraint will turn up. From 
Eq. flTO]) we can see that, in classical theory, the algebra 
of constraints is still closed. So when we write down 
the constraints with the quantum corrections, we need 
to ensure that the constraints are still closed. In other 
words, the anomaly terms, which cannot be expressed 
by the linear combination of the Hamiltonian constraint 
and the diffeomorphism constraint, should be vanished. 
On the other hand, the diffeomorphism constraint does 
not receive quantum corrections in the full theory 24] . 
so Eq.© does not change. 



The Poisson bracket between two constraints is 



{H%,D G } 



N 
1 

8^G 



(n) 
-'mil 

7 



2^ cos 

7 



(p~fk) & n (p-fk) 



d 3 xp(d c SN 3 )A 



(n)c 



(22) 
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where the anomaly part is 



7 



7 



» 



» 



^cos( pjk)& n (p-fk) - -2^ 
7 7 



(23) 



To get this, we need the Poisson bracket 

{5KZ(x),d d 6E d k (y)} = 87rG8{S d c d d 8(x,y), (24) 
{5Ei(x),d d 5K£(y)} = -^G5{5 d c d d 5{x lV ). (25) 



To cancel the anomaly part, it must be requested that 



A 



(n)c 



0, i.e. 



» 



^ = (/3+l)^cos(£ 7 fc)<S„(/nfe) 



1 - 2/3 



7 A; 



(26) 



From [23fl we know that the big bounce occur when 
: = ^, so the 
Eq.([2"6|). we have 



pc = ~, so the maximum of c^Jl is ~. Accordin; 



(/3 + l)4 Tl) /icos (Pjk) <5 n (/2 7 fc) 



J I jk 



p^k < 



(27) 



Eq. ([27|) can be seen as a limit to the evolution of cf 1 ' 
and we can restrict the range of f3 through this limit. 



Ir 



tive of the background variables 



» 



p = 2p- 



(29) 



C. Tensor mode 

In classical perturbation theory of cosmology, there is 
only one equation in tensor mode, i.e. the gravitational 
waves equation. From this equation, we know that the 
gravitational waves are massless. However, when quan- 
tum corrections are taken into account, a mass term will 
be appeared in this equation p}. It is only the conclu- 
sion calculated in the first order correction. We extend 
this method to the higher holonomy corrections, and take 
limit of n — > oo. In this way, we can find that the mass 
of gravitational waves is the intimately results from the 
quantum nature of Riemannian geometry of LQG. 

Inserting the higher holonomy corrections to Eg. dill) , 
the effective gravity part of perturbed Hamiltonian con- 
straint can be expressed as 



H%[N] 



16ttG 



d 6 xN 



2p 3 




[-&VP 



\p (SE]SE d S k c V d ) + {6Kl6K$S° k 8$ 



2 c {n l S cd S^6 e fd e E^d f E d 



VP 7 



5 3/2 



(28) 



From this Hamiltonian, one can obtain the time deriva- 



1 / c 



(«)' 



7 



cj w) P dp 
7 p dp 



c («) 

cos (p'yk) <5 n (pjk) k — 



(30) 



and the time derivative of the perturbed variable 8Ef 

SEf = [SEf, H% } = -pS a k SfSK k d - -p^hl (31 
On the other hand, one can also obtain SEf from SEf = 

i /., 



SEf = --(pht+phl 



1 ( c (n) • " 

-- [ 2p±- cos (jir/k) 0„ (Pjk) h a t + phf 



(32) 



From Eqs.([3I]) and ([32]), we have 



= - 

a 2 



/ («) (») \ 

K+h^- cos (p 7 k) e n (p 7 k) --^\ti a 



(33) 



5 



So, the 5K l a will be 
1 



8Ki = 



( >) » 
fc* + fc'ft 2-*- cos (ji 7 fc) <9„ (£ 7 fc) 



tion 

ski 



/ (n) («) \ 

2-&- cos (/2 7 /c) 0„ (/2 7 fc) 2S& J 



7 



(34) 



{SK,H matter } + 4 ( -H ^ " 5 + 2 



» / » 



» 



---^l2-^cos(/} 7 fc)6„(A7fc)-^K- (35) 



Again, we can also obtain 8K\ from Hamiltonian equa- 



From Eas.(|34p and (|35|) . one can obtain the gravitational waves equation 

» x 



fci + cos (p 7 fc) «n (A7*) ^ - o + T 



= 87rGTlQ a , 



where IIq q is the source terms from the matter Hamiltonian and 



(36) 



T, 



(n) 



Q 



2p { , 



(ri) \ 4 / (n) \ 3 2 l~ 

1 c\ ' \ COS 1 ■ 



» 



M7*0 



Q3„ (jfryfc) 



(n) 



— — c — — c 

cos (p.jk) <3„ (/i 7 fc) - J22 ^ L — cos (m/5 7 fc) <3 n (mprfk) 



1 / c 



.(«) 



» 



j |2^-7/i [sin (/2 7 fc) ©„ (£ 7 fc) - cos 2 (jx^k) <B„ (£ 7 *)] 

(c {n) - c {n) \ 2 
-1 + cos (mfr'yk) <3„ (mfljk) } - I cos (/2 7 fc) — I , 



(37) 



where 

®4A7fc)=E|P[ sin (^)] 2i " 1 . (38) 

When n = 0, 03 „ (jv)k) = 0. So it will not appear in 
conventional way. The definition of the effective mass is 



When fj/yk — > \, n — > oo, the mass term will never van- 
ish. From this, we can confirm that the nonzero mass 
of gravitational wave results from the quantum nature of 
Riemannian geometry of LQG. 

From the definition of gravitational wave [l[ , we should 

(n) 

require Tq > 0, and this is another condition to restrict 
the range of (3. 

IV. LATTICE REFINEMENTS 

In process of obtaining the range of j3 in [l| , the author 
expand the "sin" (and "cos") of Eq. (|2"6"|) fn = 0) and just 



take the first several terms of it. So in his paper, the 
first non-zero term of anomaly part is k . By this way, 
one can obtain the relationship between m and /3 from 
Eq. (|26|) fn = 0), and obtain the range of j3 by requiring 
the r£° > from Eq.([37])(n = 0). 

However, there are some problems in this method. 
First of all, the relation of m 2 = 5 + 2(3 in [l[ is obtain 
by requiring the term of k 4 is vanished. But it can not 
ensure the whole anomaly part can be canceled, because 
there are still k 6 ,k 8 , etc. Secondly, when we consider the 
higher corrections, there will be some high order terms 
like sin 3 (:r) appear. If we keep more terms of "sin" (and 
"cos"), we will find that the first non-zero term is not A: 4 , 
maybe, it will be fc 5 , it depend on how many terms you 
kept. It decides the different relation between m and /3. 
So if we want to obtain the range of /3 more accurate, we 
should not expand the "sin" (and "cos") in the equations, 
in other words, we keep all the terms of it. 

From the discussions above we can know that, Eq. ([2T|) 
and T^y* > are two restricts to /3, so we analyze these 
two restricts respectively. At first, we note that, the 
product /5 7 fc always appear together in the expression of 



6 



c mhi^ > so we set x = M7^' an ^ then is the function 

of x. When Big Bounce occurs at x = ?, and x — » with 
the expansion of the universe, we can draw the graphs of 
this function between to x = ^ with different value of 
n. 

From Fig Q] we can see that, when n = (it correspond 
the conventional holonomy corrections), the case of /3 = 



— | (which was lower bound in [l|) can fill the condition 

in) 

of c^^/i < \. But it can be more lower than that because 
(3 = —2.7 can also fill the condition. However, we will be 
concerned about the higher corrections, so let us analyze 
the case of large n. From FigfTJ we can see that, the 
larger n lead to the bigger lower bound of j3. 
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— p=-l/10 / 
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FIG. 1: The evolution of yi = -p- and y2 = c^lfi with x = [i-yk are shown in [(a),(b),(c)] and [(d),(e),(f)] respectively. The 
dash lines in (d),(e)and (f) are j/2 = §■ 



Because the Eq. (|27)l should be kept everywhere, it in- 
cludes the point of the Big Bounce i.e. x — Inserting 
x = ^ to ([2T| . it will become 



1-2/3 1 



7T 7T 

2/2-2' 



which leads to 



/?>-!■ 



(40) 



(41) 



On the other hand, the Eq. ([26|) is the relationship of 



m and /?. We insert the into Eq. (l37l) . On the right 
hand of Eq. p7p , there are still some terms which contain 
m, so we can use the following relation to replace this 



terms 



cos 



„(") 

C mh 



(n) 
C mh 

1 



(42) 

T («) 

From Eqs.(HH|), (J3Z1) and (g^l) we can see that is 
also the function of x = fijk, and there is two parameters 
n and /3 in it. 

T («) 

From the evolution of -p- with n = (see FigQJa)), 

we can see that, if we require the -p- > 0, the (3 should 
-i. With Eq.gT}, the range of f3 is 
This result is smaller than —0.1319 > 



be smaller than 

-§>/*> -I- 
/3 > — 5/2. It is because we use the "sin", not the first 
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orders of expanding term of "sin" . 

When n > 0, we find that the upper bound of j3 is 
larger than — i, and when n — > oo, the upper bound will 
be zero. We display the n = 100 and n = 1000 in Figfl] 
also. 

So, the final range of j3 should be [—1,0]. From this 
we can see that f3 — is not eliminated like in |l| and 
(3 = — | is also in this range. 

V. DISCUSSION 

In this paper, we apply the higher order holonomy cor- 
rections to the perturbation theory of cosmology. When 
we take the limit of n — > oo, the form of the LQC will be 
back to classical theory, but the effect of the quantum ge- 
ometry will be kept. From the analyses above, we know 
that the mass of gravitational wave will not be vanished 
when n — > oo . It will decrease to zero with the expansion 
of the universe. So it is the "pure" quantum effect that 
the gravitational wave have nonzero mass. 

Other important effects are related to the discrete 
space-time geometry. Discrete space means the existence 
of the area gap, and there is a function jl(p) related to 
this area gap. The form of function flip) has an impor- 
tant impact on LQC. But now the framework of theory 
is not perfect to decide this function, so we only can re- 
strict the form of the function as jj, cx pr , from some other 
aspects like effective theory and perturbation theory of 
cosmology. 

In the effective LQC framework, we apply two condi- 
tions to limit the range of f3. One is anomaly free, which 
means that the constraint algebra of vector mode should 
be closed, when we consider the the quantum effect. It is 



the mathematical requirements of the theory. This can 
restrict (3 to be [—1, +oo). 

The other condition is the requirements of positive def- 
inite mass of gravitational waves. This is the physical 
requirement. We can not ensure that the mass of grav- 
itational wave is positive when f3 > 0. And from FigfT] 
we can see that the behavior of mass between (3 > and 
(3 < is very different. Therefore it can restrict the range 
of f3 to be (— oo, 0]. This requirement seems very natural. 
However, we do not yet understand the true meaning of 
the mass of gravitational waves, so this condition is only 
an assumption. The correctness of this assumption needs 
to be verified in future studies. 

In conclusion, the range of (3 should be [—1,0]. But 
this range is only decided by perturbation theory of cos- 
mology. It cannot exclude (3 — 0. So the excluding of 
f3 = is based on the prediction of theory rather than 
theory itself. This may not be its final scope because only 
two conditions were discussed in this article. Certainly 
there are many other conditions to limit the range of pa- 
rameter. If we can restrict /3 to a unique value, say —1/2, 
from the theory itself rather than predictive power of the- 
ory, then the theory will be more self-consistency. So, in 
the future studies, we can compare different conditions 
on the parameter values to examine the self-consistency 
of theory. 
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